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Thomson scattering in non-ideal (collision-dominated) two-component plasmas is calculated ac- 
counting for electron-ion collisions as well as electron-electron correlations. This is achieved by using 
a novel interpolation scheme for the electron-electron response function generalizing the traditional 
Mermin approach. Also, ions are treated as randomly distributed inert scattering centers. The 
collision frequency is taken as a dynamic and complex quantity and is calculated from a microscopic 
~< ' quantum-statistical approach. Implications due to different approximations for the electron-electron 

correlation, i.e. different forms of the OCP local field correction, are discussed. 

G ' 
i— s 

ON 

I. INTRODUCTION 

' ' Recently, Thomson scattering has been established as a diagnostic tool for high energy laser-matter interaction in 

particular for warm dense matter [l], 0, S 3 ■ The Thomson signal probes the dynamic structure factor of the plasma 
|5(. Reversing the argument, we can synthesize the Thomson signal by using an appropriate expression for the dynamic 
structure factor and infer density and temperature conditions by matching the synthesized signal to the experimental 
one. In particular, at high densities, collisions and correlations have to be accounted for in modeling S(k,uj) Q. The 
Mermin approximation pj has been found to be a simple way of interpolating between the collision-less plasma (RPA) 
Oh. at large wavevectors and collisions in the long wavelength limit, i.e. a Drudc-likc expression for the dielectric function, 
r/} ■ However, in non-ideal plasmas, correlations beyond the RPA exist even at finite values of k and have to be taken into 
account. The traditional Mermin approach fails to incorporate these correlations. An interpolation scheme between 
static local field corrections and the Drude model by a generalized Mermin approach has been suggested in Ref. Q 
£>V using the Zubarev approach to the non-equilibrium statistical operator [§, [l(| ■ This scheme guarantees the correct 
account of electron correlations in the static limit. 

Having an interpolation scheme at our disposal, a systematic study of the influence of both electron-ion collisions 
as well as electron-electron correlations is possible. It is the objective of this paper to contribute to such a study. In 
particular, since the correct form of the dynamic as well as the static local field corrections for the interacting electron 
gas is still a matter of debate, we compare a few recent suggestions in their consequences for the Thomson scattering 
signal. 

The paper is organized as follows. In Sec. UH we give a brief review of the formalism and the approach to the 
^y-j dynamic collision frequency. Sec. lIIIl in some detail explains the different models for the dynamic local field correction 
considered here. Results for the plasmon dispersion are discussed in Sec. IIVI Finally, conclusions and an outlook 
complete this paper. 
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II. THEORETICAL BACKGROUND 
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■ We consider a neutral plasma of electrons and ions in thermal equilibrium with electron density n e , ion density 
rii = n e and temperature T. For later use, we introduce the Fermi wave vector kp = (37r 2 n e ) 1 / 3 , the Fermi energy 
Ef = h 2 k F /(2m e ), and the Brueckner parameter r s given by (47r/3)n e a|r 3 = 1, where ae is Bohr's radius. These 
parameters are relevant in our context, because we use the model of an electron gas at T = interacting with an 
inert background of ions in carrying out our exploratory calculations. 
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A. Thomson scattering and Born-Mermin approach 

It is well-known, see [1, [H, Hj] , that the experimental Thomson scattering cross section is related to the dynamic 
structure factor of all electrons in the plasma according to 



ai lau) fco 



(1) 



In this expression, cry = 6.65 x 10 - 29 m 2 is the Thomson cross section, and fco and k\ are the wavenumbers of the 
incident and the scattered light. The energy and momentum transfer are given by Ai? = Tiuj = Tuo\ — Hloo and 
Tik = Kki — Tik.0. The momentum is related to the scattering angle (L in the limit Hlo <C Hloq by fc = 47rsin(0 s /2)/Ao 
for an incident wavelength Ao . Here, we follow Chihara' s approach [5j , in that the total dynamic structure factor can 
be written in terms of contributions from free electrons , weakly and tightly bound electrons, and core electrons. In 
this paper, only the dynamic structure factor of free electrons is considered. 

In thermodynamic equilibrium, the dynamic structure factor S ee {k,Lo) and the longitudinal response function 
Xee(k,ui) are related via the fluctuation-dissipation theorem 



S ee (k,u)) = — 



1 



I _ e -huj/(k B T 



Im X e e(fc, to) 



(2) 



Theoretical approaches to dynamic structure factor of two-component plas mas have been developed starting from 
different approaches such as perturbation theory, the viscoelastic model the recurrence relation method [12j, or 
the moment approach, see Ref. 13 1. As an example for a perturbative treatment, we mention Ref. [Til. [l5|. There, 



based on the generalized linear response theory of Zubarcv, a systematic account of correlations as well as collisions 
has been accomplished by partial summation of diagram classes using thermodynamic Green's functions. While a 
detailed evaluation of the resulting expressions for the response functions is cumbersome at arbitrary wave vectors 
fc, numerical calculations have been carried out in the long-wavelength limit k — ► 0. In particular, approximative 
expressions for the collision frequency v{ui) have been studied taking care of strong collisions as well as dynamical 
screening in a consistent manner. 

To generate approximative results for the response function at finite wave vectors fc, we follow an idea suggested 
by Mermin Q. Ensuring particle number conservation by introducing local thermal equilibrium together with a 
relaxation time ansatz, the electron-electron response function is approximated by 



X {M \k,w) = Il- 



ia; \ xSf (k, lo + irj) Xcc (k, 0) 



(3) 



with a relaxation parameter ?]. For details, sec Ref. [Tof . Here, x^ PA is the electron response function in random 
phase approximation, i.e. 



x* PA (fc,-0 = 



X {0) (k,oj) 



1 - F(fc)x(°)(fc,w) 



(4) 



where x^(k,uo) is the ideal, i.c non-interacting response function , see [l6j . For T = 0, this ideal response can be 
found as 



xl 



V{k)x W ik,Lo) = -gL [g(u + z)-g(u-z 



(5) 



with u = mto /(fikkp) , z = fc/(2fc^),Xo = {^kpa-o) 1 and 



q(x) = x H — (l — x 2 ) In— — — , 
yy ' 2 1 7 x- 1 ' 



(6) 
(7) 



The function g(x) given here is a generalization of the function g(x) given in |16j to complex arguments. Furthermore 
V(k) denotes the Coulomb-potential in momentum space. Note, that Eq. ([3]) reduces to the RPA expression in the 
absence of collisions, i.c. r\ = 0. Also, in the long-wavelength limit, Eq. Q turns into the familiar Debye form, 
allowing to identify the relaxation parameter r\ as the collision frequency v{lo). 
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FIG. 1: Real part of the collision frequency u(u>) as a function of the frequency w. Various values of the Brueckner parameter 
r a are considered. 

B. Extended Mermin approach 

To account for correlations among the electrons, an extension of the traditional Mermin expression has been 
suggested, see Ref. Q. For an adiabatic model with inert ions, it reduces to replacing the RPA response function in 
Eq. ([3]) by the response function of the interacting one-component (OCP) electron gas x?<F P (k, w), 



(xm; 



(k,w) = 1- 



(8) 



where the label xM indicates the extended Mermin expression for the response function. Note, that the same expression 
has been derived independently by Barriga-Carrasco [13] ■ Traditionally, the OCP response function is represented 
using a dynamic local field correction G e e(k,u}) as 



l-l/(fc)(l-G ee (fc,u;))xi 0) (fc, W ) 



(9) 



Having a collision-less plasma v = 0, the response function is solely the OCP expression. Due to the fact, that electron- 
electron collisions do not contribute in the long- wavelength limit, i.e. G cc (k,u>) cx k 2 for k — ► 0, this expression still 
reduces to the Drudc-like form for small k with the same v(u>) as before. However, contrary to Eq. ([3]), the static 
limit is now given by the electron-electron correlation in the OCP, 



lim x 

ui—>0 



(xM) 



,OCP 



(fc.O) 



(10) 



C. Collision frequency at arbitrary degeneracy 

For the exploratory calculation discussed here, we use the collision frequency in Born approximation and for arbitrary 
degeneracy, see Ref. [l5| . 

q2 poo i 

Rc^H - - £ 2 °^ / dqq 6 V TF (q) 2 S^q) -Ime RPA , e (q,u) , (11) 

where Vtf(<z) = Viq)/^^^^, 0) is the static screened potential, Si(q) is the static structure factor of the ions, taken 
e.g. in HNC approximation or from MD simulations, f^o is a normalization volume, and CRpA.efe is the dielectric 
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FIG. 2: Static local field correction G Ee (k) for r a = 2. Approximations: Hubbard [II], Born [2§], STLS [2^], Utsumi and 
Ichimaru QJ], Farid et al. 0, Moroni et al. [H]. 

function of the electron OCP for arbitrary degeneracy as given e.g. in the article of Arista and Brandt [l6||. Again, 
we determine the RPA dielectric function for T = by using Eq. (J5j) . 

We restrict ourselves to this easily accessible expression since we want to focus on the interplay between collisions 
and electron-electron correlations and the role of different approximations for the OCP local field corrections. More 
advanced expressions are given in Ref. [l5[ and should be used for realistic calculations. Also, realistic calculations 
should include electron-electron effects on the collision frequency, which can be taken into account by increasing the 
number of moments in linear response approach, see Ref. UM as well. 

In Fig. [TJ we show the collision frequency v{ui) as a function of the frequency w for three different values of the 
Brueckncr parameter r s = 1,2,5. Note, that v(uj) is a complex quantity and the imaginary part is connected to 
Ref(w) by a Kramcrs-Kronig relation 



The account of this imaginary part is essential for obeying both, the f-sum sum rule and the perfect screening sum 
rule. While loosely speaking, the real part leads to a broadening of the plasmon at k — 0, the imaginary part produces 
a shift of the plasmon. For a static frequency uj = 0, the imaginary part vanishes, i.e. replacing the dynamic by a 
static collision frequency v(0) one ignores the shift of the plasmon position. 

III. LOCAL FIELD CORRECTIONS FOR AN INTERACTING ELECTRON GAS AT T = 



In an often used approximation, the dynamics in the local field correction is ignored reducing it to the static limit 
only, 




(12) 



A. 



Static local field correction for the OCP 





(13) 



l-V(fc)(l-Gee(fc))xHft,w) 



For the static local field correction G ee (k), a plethora of approximations have been suggested beginning with the 
original paper by Hubbard [l8||. Here, it is impossible to give an exhaustive review. Instead, we concentrate to the 
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widely- used expression of Ichimaru and Utsumi [l9| and its extension by Farid et al. [20| . For an overview of other 
approximations, see Ref. [2l|. For the sake of illustration, we compare a few approximations in Fig. H\ for r s = 2, see 
[22, [23[. Also included are results obtained with a Monte Carlo (MC) simulation by Moroni et al. [24|. Qualitatively 
similar results are obtained for other values of r s . For small values of fc, the Utsumi-Ichimaru approximation and the 
extended model of Farid et al. are identical by construction and adapted to the compressibility sum rule. There is a 
good agreement with the MC data and both, the Utsumi-Ichimaru and the Farid et al. description, while the other 
approximations do not describe these so well. In the deep inelastic regime at large fc, Farid et al. take account of the 
results by Holas [25|, that G ee (k) scales as k 2 , which is not included into the Utsumi-Ichimaru ansatz. The MC data 
seem to support this k 2 scaling. For our discussion, this difference is rather unimportant, since the plasmon ceases to 
be a well-defined mode at k Q < kp, while the differences between Ichimaru-Utsumi and Farid et al. arise for k > 2kp. 



B. Dynamic local field corrections for the OCP 



As already mentioned, there are several approximative approaches to the dynamic structure factor of the electron 
gas. Here, we use the approach of Dabrowski [26[ and the approach of Hong and Lee [13] ■ Both approaches are 
interpolation schemes incorporating sum rules and other exact properties. In particular, the static properties are 
inputs into these schemes and we can use the local field correction of Farid et al. again in this case. This would have 
been impossible if we choose the perturbative results of e.g. Richardson and Ashcroft (28j . 

Being a dynamical quantity, the local field correction is also a complex quantity obeying Kramers-Kronig like 
relations. For the real part ReG eo (fc, w) of the dynamic local field correction, the static limit is just approximated by 
G cc (k) given above. The high frequency asymptotics is given by three-frequency sum rule and therefore by the static 
structure factor and the correlated kinetic energy of the electron gas, see Ref. [2^, |3(| , 



lim G ee (k,Lu) = I(k) - 

UJ — >oo 

which in turn is given by the static structure factor S(k) 

1 



2k 2 



((Ekin) — (-Ekin)d) 



I(k) 



N 



K(k,q) (k-q) S(\q-k\)-l 



(14) 



(15) 



and 



W-k\ 2 



(16) 



Furthermore, (-Ekin) is the kinetic energy of the interacting electron gas, (-Ekin)o its non-interacting counterpart, uj p i 
is the plasma frequency. 

Also, Dabrowski incorporates the perturbative result of Glick and Long [3l| for the high-frequency behavior of the 
imaginary part of the dynamic structure factor. He extends a Pade approximation suggested by Gross and Kohn [32j 
to finite values of the wave vector, 



Im G ee (/c, to) 



a(k)uj 



(l + 6(fcV 2 ) 5 /4 



with 



»(*) 



Ck 2 



ReG oc (fc,0) — ReG ee (fc,oo) 



CDk 2 



5/3 



b(k) = 



ReG oo (fc,0) — Re G ee {k, 00) 
CDk 2 



4/3 



(17) 

(18) 
(19) 



andG = 23/66 ar si D = T(3/4)/(0i : r(5/4)) w 0.763, a = (4/(9tt)) 1 / 3 . The corresponding real part is then obtained 
by a Kramers-Kronig relation, 



/ did 
— 
-00 T 



dui' lm.G ee (k,u)') 



UJ' — LO 



(20) 
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FIG. 3: Plasmon dispersion relation in RPA approximation. Brueckner parameters r s — 1, 2, 5 are considered. The intersection 
of the plasmon dispersion with the single-particle ridge defines the wave vector ko(r s ). 



where P J indicates Cauchy principal value integration. 

As a second option for the dynamic local field correction of the OCP we introduce the above mentioned approach of 
Hong and Lee [2?l l29j|. which is based on the recurrence relation technique. Specifically, we use the lowest dynamical 
extension of the local field correction, which can be introduced by this technique. Adapting the notation to this paper, 
the dynamical local field correction reads 



G cc (k,z) = G cc (fc,0) + [G cc (fc,Oc) - Gee(M)] 



with a function cj] given by [32 



AO 



x (0) (M) 



- 11 + Af 



(21) 



(22) 



Here, the quantities A® and A° are the ideal recurrants A® = — 



A - 

L\ 2 — 



X 5 (fcy) + (fcy) 



Both approaches have been adapted in this paper to the most recent results for the static electron-electron structure 
factor S(k) and the correlated kinetic energy, see Ref. The details will be covered in an forthcoming publication 



IV. COMPARISON OF DIFFERENT APPROXIMATIONS 



Instead of showing the dynamical structure factor S(k,oj) for each of the different approximations and for various 
values of the wave vector k and the frequency tu, we introduce the plasmon position u{k) = max5 ee (fc,w) and the 

plasmon width T(q) as a signature for the influence of different effects, where we measure the width as full width 
at half maximum (FWHM). In a collision-less plasma at T = 0, the plasmon is a well-defined mode for k < ko(r s ) 
corresponding to a <5-like spike in the dynamic structure factor. For larger k, the plasmon ceases to be a well- 
defined mode and shows a broadening even in RPA, see Ref. [3. We illustrate these dispersion relations in Fig. [3] 
The wave vector k (r s ), where the plasmon dispersion V {k)x^ {k , u>) = 1 intersects with the single-particle ridge 
oj = hk (k — 2kp) /2m is smaller than /cp for all of the values r s considered here. Collisions as well as correlations 
modify these RPA dispersion relations leading to a shift and a broadening of the plasmon even at k < ko(r s ). For the 
conditions of warm dense matter, the plasmon dispersion in the traditional Born-Mcrmin approach has been studied 
by Thiclc et al. [|. 
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FIG. 4: Plasmon shift ui(k) and plasmon width T(k) as a function of the wave vector k. r a — 1 is studied. We compare the 
static local field correction to the dynamic local field correction. No collisions are included. 




0.2 0.4 0.6 0.8 1 

wave vector k [units of k p ] 



FIG. 5: Plasmon shift co(k) and plasmon width P(fc) as a function of the wave vector k. r s — 5 is studied. We compare the 
static local field correction to the dynamic local field correction. No collisions are included. 



A. Plasmon dispersion without collisions 



In a first step, we discuss the effects induced only by local field corrections, i.e. for v{u>) = 0. The modifications 
of the plasmon properties are shown in Fig. 0]for r s = 1 and in Fig. Ofor r s = 5. RPA, static local field corrections 
and dynamic local fields given by the improved Dabrowski interpolation scheme are compared. In the first case, the 
deviations of the static LFC from the RPA results arc quite small. Also, the dynamic LFC is only a minor correction 
to the static expression. This is consistent with the plasmon width, which is below 1 % for almost all wave vectors 
k. Note, that the rapid increase above 0.55 kp is due to the onset of damping in the RPA expression. For r a = 5, the 
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Brueckner parameter r s 



FIG. 6: Plasmon shift u>(0) and width F(0) as a function of the Brueckner parameter r 3 . Here, we account for electron- 
ion collisions only using the traditional Born-Mermin approximation of Eq. (0 and Eq. No local field corrections are 
considered. 



noticeable deviations occur. Also, a clear influence of dynamical LFC is visible. This is also reflected in the plasmon 
damping, where a width of up to 17 % is found. 



B. Results for the extended Mermin approach 



We start the presentation of the results for the extended Born-Mermin approach by focussing on the long-wavelength 
limit k — > 0. In this limit, the Born-Mermin ansatz reduces to a Drude-type dielectric function with a frequency- 
dependent and complex collision frequency ^(w). This frequency leads to a broadening and a shift of the plasmon as 
can be seen from the imaginary part of the inverse dielectric function, which is given by 

lme(fc^0,w) = — — ^ . (23) 

(^uj 2 — — Imzy(cj) cuj + (Kev(cu)) 2 u> 2 

Approximately, for Hev(uj) >> Im v(lo), the real part is connected to a broadening, while the imaginary part induces 
a shift of the plasmon. Local field corrections do not play any role in this limit as discussed above. The shift and 
the broadening are illustrated in Fig. [6] as a function of the Brueckner parameter r s . Since the real as well as the 
imaginary part of v{io) increase with r s for the conditions considered here, the broadening T(k = 0) of the plasmon 
also increases with r s . For large r s , the broadening is about half the size of the plasma frequency lo p \. The shift of 
the plasmon also increases with r s . However, the shift is less pronounced compared to the width of the plasmon, an 
effect consistent with Imi/(w) < Rci^(oj). 

Next, we present the wave vector dependence of the plasmon shift u(k) and the plasmon width r(fc). The results 
for the shift are shown in Fig. 003 for r s = 1,2,5. Fig. [11] displays the width for r s — 5. We compare three 
different approximations, the traditional Born-Mermin (BM) given by Eq. ([3|) and Eq. (jTTJ) , the extended Born- 
Mermin approach with static local field correction (BM+sLFC) of Eq. ((8]) together with Eq. (fT3|) and finally the 
extended Born-Mermin approach with dynamic local field corrections (BM+dLFC) by the Dabrowski ansatz, i.e. 
Eq. (TIT]) and Eq. (J^HJ) ■ Also, the RPA dispersion relation is shown, with 1 - V(k) x (0) (fc, u(k)) = for k < k a and 
the position of the maximum of Im (k , u>) for k > Icq. 

The shift at k ~ shown in these figures corresponds to the values for r s = 1, 2, 5 in Fig. [Sj In this limit, the local 
field corrections do not contribute. Thus, the different approximations merge for k — > 0. As in the long-wavelength 
limit, the deviation of all approximative expressions from the RPA results are more pronounced with increasing r s . 
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FIG. 7: Plasmon shift u{k) as a function of the wave vector k. r s — 1 is studied. BM: traditional Born-Mermin approximation. 
BM + sLFC: Born-Mermin including static local field corrections. BM + dLFC: Born-Mermin including dynamic local field 
corrections. 

As for the Born-Mermin result, it shows a systematic behavior with a switch from a blue shift to a red shift at a value 
of fc close to fco- Similar results have been reported by Thiele et al. 0], where calculations for finite temperature 
conditions at moderate degeneracy are given. We refer for details to that paper and take the traditional Born-Mermin 
results as a reference point for the inclusion of local field effects. Using the extended Born-Mermin approach together 
with a static local field correction shows a drastic change in the dispersion relation for larger values of k. This is 
expected from Fig. [21 where G ee (k) shows considerable deviations from the RPA limit, i.e. G ee (k) = 0. Also, the 
larger local field factor for increasing r s leads to a more pronounced change in the dispersion relation, as can be seen 
by comparing the three values of r s . However, once we refine the approximation by allowing for dynamic local field 
corrections, these drastic changes disappear again and a dispersion close to the original Born-Mermin curve is found, 
at least for r s — 1 and r s = 2. In the case of r s = 5, noticeable differences from both, BM and BM+sLFC remain. 
This behavior can be understood by inspection of the frequency dependence of Re G(k, uj). While Re G(k, 0) increases 
with fc, ImG(fc,oo) is considerably smaller or even negative for large k, see e.g. [301 ] . As a consequence, values at 
intermediate frequencies are considerable smaller than the static value ReG(fc,0). For larger values of fc, a zero is 
even found corresponding to a RPA-like behavior at this frequency. We illustrate this fact in Fig. [11)1 showing the real 
and the imaginary part of the dynamic local field correction as a function of the frequency lu for k = kF and r s = 2. 
Note however, that the dispersion relation is somewhat misleading. A frequency scan of the dynamic structure factor 
reveals a rich structure, which can not adequately represented by a shift and a width only. 

Finally, we investigate the plasmon width at r s — 5 for the three different approximations. Again, at k = 0, the 
broadening is only due to the collision frequency For k > fco- the rapid increase of the broadening with fc is due to 
the RPA contribution. At finite k, local field effects also contribute. In particular, the imaginary part of the dynamic 
local field correction adds to the total width of the plasmon, if one compares the BM + dLFC result to the BM curve. 
However, the net effect is not simply the sum of both contributions due to the involved arithmetics of the Mermin 
expression. This also is revealed by the BM + sLFC result, which deviates noticeably from the BM although no 
additional imaginary part for G(fc) is taken into account. 

V. CONCLUSIONS 

We devise an interpolation scheme which incorporates electron-ion collisions as well as electron-electron correlations 
by an extended Mermin approach. We apply this to an interacting electron gas at T = interacting with an inert ion 
background. As inputs act a dynamical collision frequency in Born approximation and different models for dynamic 
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wave vector k [units of k p ] 

FIG. 8: Plasmon shift u{k) as a function of the wave vector k. r s — 2 is studied. BM: traditional Born-Mermin approximation. 
BM + sLFC: Born-Mermin including static local field corrections. BM + dLFC: Born-Mermin including dynamic local field 
corrections. 

local field corrections of the interacting OCP electron gas. Plasmon properties serve as a probe for the relevance of 
collisions and correlations, respectively. 

At small wave vectors, we observe a dominance of collisions. The importance of local field corrections increase 
with increasing Brueckner parameter r s . For r s = 5 it is indispensable to account for local field correlations. Drastic 
changes in the plasmon properties which occur by taking account of static local field corrections disappear to some 
extent when using models for dynamic local field corrections. The plasmon broadening shows in general a very involved 
behavior, collisions and corrections to do not simply add up. 

To apply our model to a realistic situation of Thomson scattering in cold matter, we have to refine the input 
quantities such as the collision frequency and the model for the OCP local field corrections. At this point, there is 
certainly some improvement to be done. 

Our approach can also serve as an approach to study the interplay of impurity scattering and electron-electron 
correlation in a jcllium model description of metals. Here, connection can be made to experimental results for the 
dynamic structure factor measured by inelastic x-ray scattering (36j . 
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FIG. 10: Real and imaginary part for dynamic local field correction at r a —2 for two different wave vectors as a function of the 
frequency. 




FIG. 11: Plasmon width cu(k) as a function of the wave vector k. r a = 5 is studied. The traditional Born-Mermin without local 
field corrections is compared to the extended approach with static (BM + sLFC) and dynamic local field corrections (BM + 
dLFC). The onset ko of the single-particle excitations in RPA is also shown. 



